Conductivity of disordered quantum lattice models at infinite temperature: 

Many-body localization 



O 
(N 



I 

CZ3 



I 

o 
o 



> 

in 

cn 
o 
o 



Timothy C. Berkelbach and David R. ReichmaiQ 

Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA 

We reinvestigate the behavior of the conductivity of several disordered quantum lattice models 
at infinite temperature using exact diagonalization. Contrary to the conclusion drawn in a recent 
investigation of similar quantities in identical systems, we find evidence of a localized regime for 
strong random fields. We estimate the location of the critical field for the many-body localization 
transition for the random-field XXZ spin chain, and compare our findings with recent investigations 
in related systems. 
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I. INTRODUCTION 

Anderson localization in non-interacting systems is 
a well understood physical process whereby sufhciently 
strong disorder leads to localization of eigenstates and 
hence insulating behavior in systems that would other- 
wise be conductors.^ Extensive work performed over the 
last 50 years has established in great detail the nature 
of the Anderson transition in non-interacting systems, 
while much less work has been aimed at elucidating how 
short-ranged interactions modify the simple picture of 
transport that emerges in the non-interacting situation.^ 
This is somewhat surprising, given the fact that attention 
to the issue of interactions already appears in Anderson's 
classic 1958 paperi^ 

Recently Basko, Aleiner and Altshuler (BAA) per- 
formed a detailed diagrammatic analysis demonstrat- 
ing that weak, short-ranged electron-electron interactions 
generically lead to a finite temperature metal-insulator 
transition in systems that would be localized in the ab- 
sence of interactions.^ In fact the analysis and implica- 
tions of the work of BAA go beyond consideration of 
interacting electrons, suggesting that more general quan- 
tum entities (e.g. spins, bosons) with local interactions 
may generically fail to thermalize until a threshold energy 
is reached. The notion that such a "many-body" local- 
ization (MBL) transition may occur is at odds with the 
intuition that interactions should lead to a finite dc con- 
ductivity at all finite temperatures in analogy with the 
mechanism of phonon-mediated hopping conductivity. It 
is also at odds with, for example, the analysis of Fleish- 
man and Anderson which argues that truly short-ranged 
interactions are insufHcient to induce conductivity in an 
otherwise localized system at any temperature^ 

The work of BAA has motivated more recent investi- 
gations of the possibility of a finite temperature transi- 
tion from a localized (non-ergodic) phase to a delocal- 
ized phase where interactions afford thermalization of 
the system. These works, both analytical and numeri- 
cal, have reached somewhat conflicting conclusions i^"— 
In this work, we reconsider the analysis of Karahalios et 
al.-' Via examination of the conductivity, these authors 
concluded that, in general, finite temperature systems of 



one-dimensional interacting spins are always conducting, 
thus contradicting the claim of BAA. 



II. IMAGINARY FREQUENCY 
CONDUCTIVITY 

Here, as in previous work, we make use of the impor- 
tant observation of Oganesyan and Huse that the many- 
body localization transition may be probed at infinite 
temperature by varying the disorder strength.^ This sim- 
plifies the problem by reducing the number of control pa- 
rameters that may be varied to tune the system from a 
delocalized to a localized phase. As in the work of Kara- 
halios et al., we examine one-dimensional spin chains via 
exact diagonalization, calculating the conductivity via 
the Kubo formula. An important conclusion of our work 
is that sufficient care needs to be exercised in the inter- 
pretation of the zero frequency conductivity as a function 
of disorder strength and level broadening. 

For a finite system of length L at T — >■ oo, the Kubo 
formula for the conductivity is given by 



cr(w) 



— lim 

L n^o 



e^i'^+^^)'{j{t)m)dt, (1) 



where /3 = l/ksT, j{t) is the current operator at time t, 
and 77 can be thought of as both a numerical tool for con- 
vergence as well as a phenomenological level broadening 
for discrete spectraJ^ The real part of the conductivity 
may be decomposed as cr'(a;) ~ D5{uj) + CTiog(i^), where 
the Drude weight D measures purely ballistic conduction 
and arises due to pairs of degenerate states connected 
by the current operator. However, in the systems stud- 
ied here, all level degeneracies are lifted in the presence 
of disorder, and conductivity is purely diffusive. Thus 
we may take as our definition of the dc conductivity 
CTdc = cr(a; 0) without concern for the Drude con- 
tribution. By employing the spectral representation of 
the Hamiltonian we arrive at the 'imaginary frequency' 
dc conductivity, 

^(^.) = ^ElH.NI^;^,:p^, (2) 



where H\n) = En\n), u}„ 
partition function. 



En — Em, and Z is the 



The authors of Ref.[8] calculate the full frequency- 
dependent (ac) conductivity spectrum using a level- 
broadening binning procedure and draw conclusions re- 
garding dc conductivity based on the w — > behavior. 
However, as we will show, all finite-sized systems with 
level broadening will exhibit a non-zero dc conductivity, 
and thus such an analysis is inconclusive. Rather, it is the 
conductivity's dependence on this level broadening which 
allows one to draw conclusions regarding conducting and 
insulating behaviors in the thermodynamic limit. 

As discussed by Thouless and Kirkpatrickfi^ finite sys- 
tems should display a simple asymptotic behavior for the 
dc conductivity that scales as ry for small ry and rj~^ for 
large 77. The distinction between conductor and insula- 
tor manifests in the behavior between these asymptotic 
regimes. In particular, we expect that a system with 
insulating behavior will exhibit an imaginary frequency 
dc conductivity with well resolved 77 and rj~^ regimes 
separated by a simple maximum. Finite sized systems 
expected to behave as conductors in the L — > 00 limit 
exhibit a broad crossover between these regimes, with a 
plateau signifying the onset of a true dc conductivity, ^iii^ 
Although the dc conductivity is well defined only in the 
thermodynamic limit (specifically, L — 00, then 77 — >■ 0), 
we see evidence for dc conductivity manifesting itself even 
at the small system sizes accessible by exact diagonaliza- 
tion. 

It should be pointed out that although the analysis em- 
ployed here was originally developed for non-interacting 
systems, our results empirically show that it is equally 
applicable to interacting ones, by replacing the single- 
particle levels with many-body levels. Specifically, we 
locate an insulating regime in which <j{ivi) on rj when 77 is 
less than the level spacing in the many-body localization 
volume, which does not scale with the size of the sys- 
tem. This behavior is to be contrasted with an observed 
metallic regime, which has (7(777) 77 as long as 77 is less 
than the many-body level spacing in the system volume 
- a spacing which vanishes in the thermodynamic limit 
yielding a dc conductivity plateau at small to intermedi- 
ate 1]. 



III. QUANTUM LATTICE MODELS 



We study two quantum lattice models in the presence 
of disorder: the XX Z spin chain and the t-t'-V model of 
spinless fermions, originally studied in its disordered form 
by Oganesyan and Huse^ and more recently by Monthus 
and GarelflS The disordered XXZ chain is given by the 
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FIG. 1: Local spin-spin correlation function, {Sj {t)Sj (0)) for 
the disordered XXZ chain with W = 5. Results are shown 
for the non-interacting case, A — 0.0 (filled circles) and for 
the interacting case, A = 0.5 (solid line). Error bars for the 
non-interacting case (not shown for clarity) are of the same 
order as those shown for the interacting case. 



Hamiltonian 



where we choose the random fields Wj uniformly from 
[—IF, IF]. The current operator for the XXZ chain is 
given by 

L 

jxxz^jJ2[s;;s]^,-sys^+^. (4) 

The disordered t-t'-V model is described by the Hamil- 
tonian 



Ht-t'-v - 



- 1 [CjCj+i + cy^j^Cj 

- t' {c]cj+2 + c]+2Cj) 



(5) 



+ V{n 



1 



where, following Oganesyan and Huse, the random on- 
site energies Wj are chosen from a Gaussian distribution 
with mean and variance IF^. The t-t'-V model's cur- 
rent operator is 



Jt-t'-V 



(6) 



2f 



+2 - C]+2Cj 
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r|, imaginary part of frequency [J] 

FIG. 2: The dc conductivity as a function of the imaginary 
frequency, rj, for the disordered XXZ chain (J — 1) with 
W = 1 — 6 [(a) — (f)] for the non- interacting case, A = 0.0 
(black hne, circles) and for two interacting cases, A = 0.5 
(red line, squares) and A = 1.0 (green line, diamonds). Error 
bars are smaller than the symbols. The blue line shows linear 
slope. 
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FIG. 3: System size dependence of the dc conductivity as 
a function of the imaginary frequency, rj, for the disordered 
XXZ chain (J = 1) with A = 0.5 and W = 1. Conductivities 
are presented for L = 10 (green plusses), 12 (blue diamonds), 
14 (red squares), and 16 (black circles), showing the develop- 
ment of the dc conductivity plateau in the L — > oo limit. 



spin correlation exhibits quantitatively similar relaxation 
compared to the localized system and shows no sign of de- 
cay from a plateau (the analog of the Edwards- Anderson 
parameter), indicative of a glassy phase. 



IV. RESULTS AND THE MANY-BODY 
LOCALIZATION TRANSITION 



In what follows, we restrict our Hilbert space to S^^f = 
for the XXZ chain and to half-filling for the t-t'-V 
model. We have explored other alternatives and find our 
results to be qualitatively similar. All results are pre- 
sented for L = 14 with periodic boundary conditions, 
although we have studied systems as large as L = 16 
and find our conclusions unaltered. We average over 
Nr = 100 independent realizations of disorder and cal- 
culate error bars as the standard deviation of the mean, 
^al\fWr, where a (not to be confused with the conduc- 
tivity) is the standard deviation across disorder realiza- 
tions. Our results appear to be converged, although for 
such a small range of system sizes one should view such 
statements with care. This is especially true for smaller 
values of disorder, where finite localization lengths, if 
they exist, can clearly be larger than the system sizes 
accessible from exact diagonalization. 

The aforementioned claim of Karahalios et al. is rather 
surprising in light of the fact that an earlier tDMRG 
calculation performed on strongly disordered spin chains 
found evidence of a localized phase, at least when viewed 
from the standpoint of the local spin-spin correlation 
functionii In Fig. [T] we reproduce this result (for a smaller 
system and the shorter time scales accessible in exact di- 
agonalization) and compare it with dynamics in the non- 
interacting system which is known to be localized. It 
is clear that, for the interacting system, the local spin- 



We turn next to an investigation of the conductivity in 
the above disordered lattice models. As discussed above, 
it is useful to compare the r] dependence of the dc con- 
ductivity in the interacting models directly with their 
noninteracting counterparts to set a baseline for local- 
ization. In Fig. [5] we present results for two different 
values of the anisotropy A in the XXZ spin chain. We 
extend the results of Karahalios et al., who examine dis- 
order strengths only as high as = 1, by performing 
our calculations up to = 6. Clearly, by = 4 in 
both interacting cases the conductivity curves are essen- 
tially indistinguishable from the non-interacting localized 
case. Thus, at least with regard to exact diagonalization 
on these system sizes and based solely on examination of 
the conductivity, the interacting behavior is identical to 
the non-interacting, localized behavior. Using this con- 
dition, we can place an upper limit, Wc < 4 in both 
cases A = 0.5 and A = l.Oi^i It should be noticed as 
well that already at IF = 1 there is significant structure 
in (T(?ry), exhibiting a nearly flat region in between the 
small and large r\ regimes. This suggests that these in- 
teracting data are in a conducting regime, although care 
must be used because this also might be an indication 
of localization behavior on length scales larger than we 
can access via exact diagonalization. One should also 
note that this conducting behavior is fully consistent 
with the conclusions of Karahalios et al. a.t W = 1.0. 
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FIG. 4: The same as in Fig[2] but for tlie disordered t-t'-V 
model (t^t' = 1) witii W = 3 (a), W = 5 (b), W ^ 10 (c), 
and W = 16 (d), for tiie non-interacting case, V = (black 
line, circles) and for the interacting case, V = 2 (red line, 
squares). Where not shown, error bars are smaller than the 
symbols. The blue line shows linear slope. 



We have also examined the behavior of the disordered 
t-t'-V model of Oganesyan and Huse, the results of which 
are presented in Fig. H) The behavior is qualitatively the 
same as above, suggesting the existence of a MBL transi- 
tion in this model as well. However, it should be pointed 
out that the crossover between apparent conducting and 
insulating behavior in cr^ir]) is significantly broader in the 
t-t'-V model than in the XX Z systems, making a predic- 
tion of the critical disorder strength difficult. At larger 
system sizes this crossover should become more abrupt, 
but unfortunately such sizes are beyond the reach of ex- 
act diagonalization. Despite the above difficulties, we 
would expect, based on the same means of analysis pre- 
sented above, that the critical value of W in this model 
is higher than the Wc ~ 5 range found in the real-space 
renormalization group calculation of Monthus and Garel. 
Although, it is not clear when studying finite systems 
that different quantities, such as those investigated here 
and by Monthus and Garel, should behave in a similar 
manner. However, our result does strongly suggest con- 
ducting behavior at W ^ 5. 



With the above caveats, we can place the critical value of 
W for a many-body localization transition in the range 
3 < Wc < 4ji^ It is thus unsurprising that Karahalios 
et al. found no evidence for the MBL transition, as we 
have shown that it occurs at disorder strengths larger 
than those investigated in Ref.[8]. We have additionally 
analyzed adjacent many-body level-spacings (not shown 
here) , whose crossover from the Gaussian Orthogonal En- 
semble to Poisson statisticsi^""— occurs at this same criti- 
cal strength of disorder, confirming the robustness of our 
approach. 

In order to confirm our expectations of a conducting 
phase in the thermodynamic limit, L — >■ cxi, we have ex- 
amined the effects of system size on the conductivity of 
the interacting XX Z spin chain with A = 0.5. We focus 
on the disorder strength W = 1 because of its apparent 
conducting behavior in Fig. [2] Furthermore, one is in 
danger of approaching localization lengths equal to the 
size of the sytem for values of disorder much smaller than 
this. We show in Fig.[3]the conductivities for system sizes 
L — 10, 12, 14, and 16. Results for L = 16 are shown 
for Nr = 50 realizations of disorder. Clearly, the plateau 
becomes more resolved for larger system sizes, strongly 
suggesting a non-zero dc conductivity in the thermody- 
namic limit. Furthermore, the plateau's growth extends 
toward smaller rj in agreement with the many-body level 
spacing discussion above. 



V. CONCLUSIONS 

To summarize, we have carried out a systematic in- 
vestigation of the diffusive conductivities of two common 
disordered quantum lattice models: the XX Z spin chain 
and the t-t'-V model of spinless fermions. We find that 
by studying the behavior of cr(i??) we can place reason- 
able bounds on the location of the MBL transition. By 
examining disorder strengths higher than those explored 
in previous works, wc find for the disordered XX Z chain 
(both A — 0.5 and A = 1.0) that the finite size conduc- 
tivity extrapolates quantitatively to the non-interacting 
values by « 4. Our results are also qualitatively con- 
sistent with the recent work of Monthus and Garel, but 
we would expect based on finite size studies of the con- 
ductivity, that the transition occurs at a disorder value 
larger than that found by their numerical renormaliza- 
tion group procedure. 
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It should be noted that it is possible that a residual conduc- 
tivity, impossible to resolve via an exact diagonalization 
study of (Tdc on small systems, exists. We see no evidence 
for finite size effects here, though we are restricted to very 
small L. 

After this work was completed we became aware of the 
work of Pal and Huse (arXiv:1003.2613vl), who also in- 
vestigate the XXZ chain with A = 1. Their conclusions 
are quantitatively consistent with our results. It should be 
noted that their work goes beyond that presented here by 
connecting the MBL transition to the infinite randomness 
universality class. 



